Epigenetic regulation of the honey bee transcriptome: unravelling the nature of methylated genes

Background: Epigenetic modification of DNA via methylation is one of the key inventions in eukaryotic evolution. It provides a source for the switching of gene activities, the maintenance of stable phenotypes and the integration of environmental and genomic signals. Although this process\ud is widespread among eukaryotes, both the patterns of methylation and their relevant biological roles not only vary noticeably in different lineages, but often are poorly understood. In addition, the evolutionary origins of DNA methylation in multicellular organisms remain enigmatic. Here we used a new 'epigenetic' model, the social honey bee Apis mellifera, to gain insights into the significance of methylated genes.\ud \ud Results: We combined microarray profiling of several tissues with genome-scale bioinformatics and bisulfite sequencing of selected genes to study the honey bee methylome. We find that around 35% of the annotated honey bee genes are expected to be methylated at the CpG dinucleotides by a highly conserved DNA methylation system. We show that one unifying feature of the methylated genes in this species is their broad pattern of expression and the associated 'housekeeping' roles. In contrast, genes involved in more stringently regulated spatial or temporal functions are predicted to be un-methylated.\ud \ud Conclusion: Our data suggest that honey bees use CpG methylation of intragenic regions as an epigenetic mechanism to control the levels of activity of the genes that are broadly expressed and might be needed for conserved core biological processes in virtually every type of cell. We discuss the implications of our findings for genome-scale regulatory network structures and the evolution\ud of the role(s) of DNA methylation in eukaryotes. Our findings are particularly important in the context of the emerging evidence that environmental factors can influence the epigenetic settings of some genes and lead to serious metabolic and behavioural disorders

An efficient approach to screening epigenome-wide data

Screening cytosine-phosphate-guanine dinucleotide (CpG) DNA methylation sites in association with some covariate(s) is desired due to high dimensionality. We incorporate surrogate variable analyses (SVAs) into (ordinary or robust) linear regressions and utilize training and testing samples for nested validation to screen CpG sites. SVA is to account for variations in the methylation not explained by the specified covariate(s) and adjust for confounding effects. To make it easier to users, this screening method is built into a user-friendly R package, ttScreening, with efficient algorithms implemented. Various simulations were implemented to examine the robustness and sensitivity of the method compared to the classical approaches controlling for multiple testing: the false discovery rates-based (FDR-based) and the Bonferroni-based methods. The proposed approach in general performs better and has the potential to control both types I and II errors. We applied ttScreening to 383,998 CpG sites in association with maternal smoking, one of the leading factors for cancer risk

Identifying heterogeneous transgenerational DNA methylation sites via clustering in beta regression

This paper explores the transgenerational DNA methylation pattern (DNA methylation transmitted from one generation to the next) via a clustering approach. Beta regression is employed to model the transmission pattern from parents to their offsprings at the population level. To facilitate this goal, an expectation maximization algorithm for parameter estimation along with a BIC criterion to determine the number of clusters is proposed. Applying our method to the DNA methylation data composed of 4063 CpG sites of 41 mother–father-infant triads, we identified a set of CpG sites in which DNA methylation transmission is dominated by fathers, while at a large number of CpG sites, DNA methylation is mainly maternally transmitted to the offspring

Evaluating the efficacy of breastfeeding guidelines on long-term outcomes for allergic disease

Background: WHO guidelines advocate breastfeeding for six months, and EAACI recommends exclusive breastfeeding for 4-6 months. However, evidence for breastfeeding to prevent asthma and allergic disease is conflicting. We examined whether following recommended breastfeeding guidelines alters the long-term risks of asthma, eczema, rhinitis, or atopy.Methods: The effect of non-exclusive (0, >0-6, >6 months), and exclusive breastfeeding (0, >0-4, >4 months) on repeated measures of asthma (10, 18 years), eczema, rhinitis, and atopy (1-or-2, 4, 10,18 years) risks were estimated in the IoW cohort (n=1456) using log-linear models with generalised estimating equations. The Food Allergy and Intolerance Research (FAIR) cohort (n=988), also from the IoW, was examined to replicate results.Results: Breastfeeding (any or exclusive) had no effect on asthma and allergic disease in the IoW cohort. In the FAIR cohort, any breastfeeding for >0-6 months protected against asthma at 10 years (RR=0.50, 95%CI=0.32-0.79, p=0.003) but not other outcomes, while exclusive breastfeeding for >4 months protected against repeated rhinitis (RR=0.36, 95%CI=0.18-0.71, p=0.003). Longer breastfeeding was protective against late-onset wheeze in the IoW cohort.Conclusion: The protective effects of non-exclusive and exclusive breastfeeding against long-term allergic outcomes were inconsistent between these co-located cohorts, agreeing with previous observations of heterogeneous effects. Although breastfeeding should be recommended for other health benefits, following breastfeeding guidelines did not appear to afford consistent protection against long-term asthma, eczema, rhinitis or atopy. Further research is needed into the long-term effects of breastfeeding on allergic disease

Tetanus vaccination is associated with differential DNA-methylation: Reduces the risk of asthma in adolescence

BackgroundVaccinations have been suggested to be associated with increased risk of allergic diseases. Tetanus vaccination is one of the most frequently administered vaccines as a part of wound management and was also found to be associated with increased serum IgE levels. We hypothesized that the vaccination modifies the risk of allergic diseases through epigenetic changes such as DNA methylation.MethodData on tetanus vaccination between 10 and 18 years of age was collected from a birth cohort established on the Isle of Wight UK in 1989. DNA methylation data were collected from individuals at different ages (at birth [n = 30], age 10 [n = 34], age 18 [n = 245] and during pregnancy [n = 121]) using the Illumina Infinium HumanMethylation450 K array. Firstly, we performed an epigenome-wide screening to identify cytosine-phosphate-guanine sites (CpGs) associated with tetanus vaccination in 18-year-olds. Secondly, we tested their association with asthma, allergic sensitization, eczema, serum IgE and pulmonary lung function (FVC, FEV1, FEV1/FVC, and FEF25-75%). We then described changes in the methylation of the selected CpG sites over age, and by vaccination status.ResultsTetanus vaccination was found to be associated with decreased methylation of cg14472551 (p value 0.5 × 10?5, FDR-adjusted p value 2.1 × 10?4) and increased methylation of cg01669161 (p value 0.0007, FDR-adjusted p value 0.014). Both CpGs, in turn, were associated with decreased risk of asthma at 18 years of age. Cg14472551 is located in an intron of KIAA1549L, whose protein binds to a B-cell commitment transcription factor; cg01669161 is located between an antisense regulator of the proteasome assembly chaperone PSMG3, and TFAMP1, a pseudogene. Increased methylation of cg01669161 was also associated with decreased serum IgE levels.ConclusionDNA methylation changes following tetanus vaccination may offer a novel prospect to explain a differential occurrence of asthma in adolescence

The interplay of DNA methylation over time with Th2 pathway genetic variants on asthma risk and temporal asthma transition

BackgroundGenetic effects on asthma of genes in the T-helper 2 (Th2) pathway may interact with epigenetic factors including DNA methylation. We hypothesized that interactions between genetic variants and methylation in genes in this pathway (IL4, IL4R, IL13, GATA3, and STAT6) influence asthma risk, that such influences are age-dependent, and that methylation of some CpG sites changes over time in accordance with asthma transition. We tested these hypotheses in subsamples of girls from a population-based birth cohort established on the Isle of Wight, UK, in 1989.ResultsLogistic regression models were applied to test the interaction effect of DNA methylation and SNP on asthma within each of the five genes. Bootstrapping was used to assess the models identified. From 1,361 models fitted at each age of 10 and 18 years, 8 models, including 4 CpGs and 8 SNPs, showed potential associations with asthma risk. Of the 4 CpGs, methylation of cg26937798 (IL4R) and cg23943829 (IL4) changes between ages 10 and 18 (both higher at 10; P?=?9.14?×?10?6 and 1.07?×?10?5, respectively).At age 10, the odds of asthma tended to decrease as cg12405139 (GATA3) methylation increased (log-OR?=??12.15; P?=?0.049); this effect disappeared by age 18. At age 18, methylation of cg09791102 (IL4R) was associated with higher risk of asthma among subjects with genotype GG compared to AG (P?=?0.003), increased cg26937798 methylation among subjects with rs3024685 (IL4R) genotype AA (P?=?0.003) or rs8832 (IL4R) genotype GG (P?=?0.01) was associated with a lower asthma risk; these CpGs had no effect at age 10. Increasing cg26937798 methylation over time possibly reduced the risk of positive asthma transition (asthma-free at age 10???asthma at age 18; log-OR?=??3.11; P?=?0.069) and increased the likelihood of negative transition (asthma at age 10???asthma-free at age 18; log-OR?=?3.97; P?=?0.074).ConclusionsThe interaction of DNA methylation and SNPs in Th2 pathway genes is likely to contribute to asthma risk. This effect may vary with age. Methylation of some CpGs changed over time, which may influence asthma transition

Gene expression differences in relation to age and social environment in queen and worker bumble bees

Eusocial insects provide special insights into the genetic pathways influencing aging because of their long-lived queens and flexible aging schedules. Using qRT-PCR in the primitively eusocial bumble bee Bombus terrestris (Linnaeus), we investigated expression levels of four candidate genes associated with taxonomically widespread age-related pathways (coenzyme Q biosynthesis protein 7, COQ7; DNA methyltransferase 3, Dnmt3; foraging, for; and vitellogenin, vg). In Experiment 1, we tested how expression changes with queen relative age and productivity. We found a significant age-related increase in COQ7 expression in queen ovary. In brain, all four genes showed higher expression with increasing female (queen plus worker) production, with this relationship strengthening as queen age increased, suggesting a link with the positive association of fecundity and longevity found in eusocial insect queens. In Experiment 2, we tested effects of relative age and social environment (worker removal) in foundress queens and effects of age and reproductive status in workers. In this experiment, workerless queens showed significantly higher for expression in brain, as predicted if downregulation of for is associated with the cessation of foraging by foundress queens following worker emergence. Workers showed a significant age-related increase in Dnmt3 expression in fat body, suggesting a novel association between aging and methylation in B. terrestris. Ovary activation was associated with significantly higher vg expression in fat body and, in younger workers, in brain, consistent with vitellogenin's ancestral role in regulating egg production. Overall, our findings reveal a mixture of novel and conserved features in age-related genetic pathways under primitive eusociality

Epigenome-wide association study of asthma- and wheeze characterizes loci within HK1

Background:To identify novel epigenetic markers of adolescent asthma and replicate findings in an independent cohort, then explore whether such markers are detectable at birth, predictive of early-life wheeze, and associated with gene expression in cord blood.Methods:We performed epigenome-wide screening with recursive random forest feature selection and internal validation in the IOW birth cohort. We then tested whether we could replicate these findings in the independent cohort ALSPAC and followed-up our top finding with children of the IOW cohort.Results:We identified 10 CpG sites associated with adolescent asthma at a 5% false discovery rate (IOW, n = 370), five of which exhibited evidence of associations in the replication study (ALSPAC, n = 720). One site, cg16658191, within HK1 displayed particularly strong associations after cellular heterogeneity adjustments in both cohorts (ORIOW = 0.17, 95% CI 0.04–0.57) (ORALSPAC = 0.57, 95% CI 0.38–0.87). Additionally, higher expression of HK1 (OR = 3.81, 95% CI 1.41–11.77) in cord blood was predictive of wheezing in infancy (n = 82).Conclusion:We identified novel associations between asthma and wheeze with methylation at cg16658191 and the expression of HK1, which may serve as markers of, predictors of, and potentially etiologic factors involved in asthma and early life wheeze

The genomes of two key bumblebee species with primitive eusocial organization

Background: The shift from solitary to social behavior is one of the major evolutionary transitions. Primitively eusocial bumblebees are uniquely placed to illuminate the evolution of highly eusocial insect societies. Bumblebees are also invaluable natural and agricultural pollinators, and there is widespread concern over recent population declines in some species. High-quality genomic data will inform key aspects of bumblebee biology, including susceptibility to implicated population viability threats. Results: We report the high quality draft genome sequences of Bombus terrestris and Bombus impatiens, two ecologically dominant bumblebees and widely utilized study species. Comparing these new genomes to those of the highly eusocial honeybee Apis mellifera and other Hymenoptera, we identify deeply conserved similarities, as well as novelties key to the biology of these organisms. Some honeybee genome features thought to underpin advanced eusociality are also present in bumblebees, indicating an earlier evolution in the bee lineage. Xenobiotic detoxification and immune genes are similarly depauperate in bumblebees and honeybees, and multiple categories of genes linked to social organization, including development and behavior, show high conservation. Key differences identified include a bias in bumblebee chemoreception towards gustation from olfaction, and striking differences in microRNAs, potentially responsible for gene regulation underlying social and other traits. Conclusions: These two bumblebee genomes provide a foundation for post-genomic research on these key pollinators and insect societies. Overall, gene repertoires suggest that the route to advanced eusociality in bees was mediated by many small changes in many genes and processes, and not by notable expansion or depauperation